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ABSTRACT 

We report the identification of the optical afterglow of GRB 980613 in R- and 
/-band images obtained between 16 and 48 hours after the gamma-ray burst. 
Early near-infrared (NIR) H and K' observations are also reported. The after- 
glow was optically faint (R ~ 23) at discovery but did not exhibit an unusually 
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rapid decay (power-law decay slope a < 1.8 at 2a). The optical/NIR spec- 
tral index (/3rh < 1-1) was consistent with the optical-to- X-ray spectral index 
(Aix ~ 0.6), indicating a maximal reddening of the afterglow of ~ 0.45 mag in 
R. Hence the dimness of the optical afterglow was mainly due to the fairly flat 
spectral shape rather than internal reddening in the host galaxy. We also present 
late-time HST/STIS images of the field in which GRB 980613 occurred, obtained 
799 days after the burst. These images show that GRB 980613 was located close 
to a very compact, blue V = 26.1 object inside a complex region consisting of 
star-forming knots and/or interacting galaxy fragments. Therefore, GRB 980613 
constitutes a strong case for the association of cosmological gamma-ray bursts 
with star-forming regions. 

Subject headings: cosmology: observations — gamma rays: bursts — galaxies: 
star formation — galaxies: photometry 



1. Introduction 

Gamma-ray bursts (GRBs) are the most energetic events in the Universe. The currently 
favored models for long-duration GRBs (lasting more than about one second) involve the 
collapse of a very massive progenitor star (see e.g., MacFadyen & Woosley (1999), Vietri & 
Stella (1998), and Bloom, Kulkarni & Djorgovski (2002a)). Because of the short-lived nature 
of massive stars, GRBs are believed to originate in regions where stars are born. Hence, in 
this scenario, GRBs are expected to be efficient tracers of massive-star formation in galaxies 
(Paczyhski 1998). Consequently, GRBs may trace the (massive) star- formation rate in the 
Universe as a function of redshift (Totani 1997; Wijers et al. 1998; Blain & Natarajan 2000). 
Evidence supporting this hypothesis includes (i) the association of SN1998bw and GRB 
980425 which occurred in a star-forming region (Galama et al. 1998; Fynbo et al. 2000), (ii) 
the association of cosmological GRBs with supernovae (e.g., Bloom et al. 1999; Galama et 
al. 2000; Bjornsson et al. 2001; Bloom et al. 2002; Garnavich et al. 2002), (iii) the discovery 
of X-ray absorption edges and lines in GRB afterglows (Amati et al. 2000; Piro et al. 2000; 
Reeves et al. 2002), (iv) the locations of GRBs close to star- forming regions in interacting 



1 Bascd on observations with the Nordic Optical Telescope, which is operated on the island of La Palma 
jointly by Denmark, Finland, Iceland, Norway and Sweden, at the Spanish Observatorio del Rocque de los 
Muchachos of the Instituto de Astrofisica de Canarias, and on observations with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc. under NASA contract NAS5-26555. 
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galaxies (Holland & Hjorth 1999), in starburst or ultraluminous infrared galaxies (ULIRGs) 
(Djorgovski et al. 1998a; Holland et al. 2001; Berger, Kulkarni & Frail 2001) or in Lya- 
emitting galaxies (Fynbo et al. 2002), (v) indications of circumburst density distributions 
similar to those generated by stellar winds (Jaunsen et al. 2001), (vi) the existence of GRBs 
at very high redshifts (Andersen et al. 2000), and (vii) the inferred high column densities of 
neutral Hydrogen along the line of sight to GRBs (Wijers & Galama 2001). 

GRB 980613 triggered the BeppoSAX Gamma-Ray Burst Monitor (40-700 keV, Costa 
et al. (1998a)) on 1998 June 13 at 04:51:06 UT and was simultaneously detected and posi- 
tioned at the equatorial coordinates R.A. (J2000.0) = 10hl7m44s, Dec. (J2000.0) = +71°29'9 
(4' error radius) by the Wide Field Cameras instrument (2-26 keV, Jager et al. (1997)) on 
board the same satellite (Smith et al. 1998). The event was also detected by BATSE/ CGRO 
with a flux and fluence near the middle of the BATSE burst distribution (Woods, Kippen 
& Connaughton 1998). As measured by BeppoSAX, it was composed of a single pulse 
of duration « 15 s in 7-rays and « 50 s in X-rays. The peak flux and fluence were 
(1.3 ± 0.4) x 10~ 7 erg cm" 2 s" 1 and (1.0 ± 0.2) x 10~ 6 erg cm" 2 in the 40-700 keV energy 
band, and (2.0 ± 0.3) x 10~ 8 erg cm" 2 s" 1 and (2.7 ± 0.4) x 10~ 7 erg cm" 2 in the 2-26 keV 
energy band. The BeppoSAX Narrow Field Instruments (Boella et al. 1997) were pointed 
to the GRB 980613 error box about 9 h after the burst for a total time of about 100 ks, 
leading to the discovery of a possible X-ray afterglow source at R.A. (J2000.0) = 10hl7m53s, 
Dec. (J2000.0) = +71°27'24" (1' error radius) with an average flux of (1.1 ± 0.3) x lO" 13 
erg cm" 2 s" 1 (Costa et al. 1998b). GRB 980613 is the faintest burst detected by BeppoSAX 
in the 7-ray band but is relatively rich in X-rays (Costa et al. 1998c; Costa 1999). It belongs 
to the less luminous X-ray class of Boer & Gendre (2000). 

Following the discovery of the optical afterglow by Hjorth et al. (1998) (described in 
more detail in this paper) an object identified as the host galaxy was detected by Djorgovski 
et al. (1998b) and suggested to be part of an interacting system by Hjorth et al. (1999). The 
complex region in which GRB 980613 occurred was studied in more detail by Djorgovski et al. 
(2000) who concluded that the burst was located in a starburst galaxy at z = 1.0969±0.0002 
which is interacting with one or more of several galaxies or galaxy fragments in the vicinity 
of the GRB. This led Djorgovski et al. (2000) to propose that the GRB originated from a 
merger-induced starburst. 

We here present the initial identification and optical/near- infrared follow-up observa- 
tions of the afterglow which are used to constrain its exact location, decay rate and spectral 
index (§2). We next discuss the properties of the host and its environment, as observed with 
the Hubble Space Telescope (HST) (§3). We conclude by discussing the implications for the 
nature of dark bursts and for the birth places of GRBs. 
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We assume a cosmology where H = 70 km s _1 Mpc" 1 , Vl m = 0.3, and Qa = 0.7. In 
this case a redshift of 1.0969 corresponds to a luminosity distance of 7.41 Gpc and a distance 
modulus of 44.35. One arcsecond corresponds to 8.2 proper kpc, and the lookback time is 
8.1 Gyr. 

2. The optical afterglow 

The optical and near-infrared observations are reported in Table 1. The discovery of the 
optical afterglow was made on the basis of a comparison of the 1998 June 13.9 UT and 17.9 
UT NOT images as described by Hjorth et al. (1998) and shown in Fig. 1. There were only 
two detections of the optical afterglow in each of the R and I bands from NOT and MDM. In 
the near infrared the afterglow was not detected to a limit of H = 20.9, but a faint object was 
however detected with K' = 18.8 ± 0.3. The magnitude of the host galaxy is K = 19.1 ± 0.2 
(R. Chary, private communication, see also Chary, Becklin & Armus (2002)). Even taking 
the offset between K and K' into account this indicates a substantial contribution from the 
host galaxy. For this reason we disregard the K' data in the remainder of this paper. 

Djorgovski et al. (2000) reported the detection of several objects within a few arcseconds 
of the afterglow. Using the HST data presented in §3 we subtracted the contribution from 
these objects to derive the corrected optical afterglow magnitudes presented in Table 1. 
At discovery, 16 hours after the burst, the brightness of the optical afterglow was R = 
23.19 ± 0.20, making it one of the faintest known optical afterglows (Fynbo et al. 2001; 
Lazzati, Covino & Ghissellini 2002; Bloom et al. 2002c). 

An astrometric fit to the Guide Star Catalog II yielded the afterglow coordinates 
R.A. (J2000.0) = 10hl7m57.90s, Dec. (J2000.0) = +71°27'26 / . , 3 with an rms uncertainty 
of ±0"27. 

As observed for most afterglows (and expected from theory cf. e.g. Meszaros, Rees & 
Papathanassiou (1994); Sari, Piran & Narayan (1998)), the optical spectrum and light decay 
of the optical afterglow can be described by a power-law, i.e. f v oc (t — ^grb) - a ^ ^- Formally, 
the decay slope derived from the i?-band data (16-24 hours after the burst) is an = 0.8 ±0.5 
while ai = 1.3 ± 0.8 from the J-band data (24-48 hours after the burst). These values 
are consistent with the decay slope of 1.05 ± 0.37 observed in the X-rays (Boer & Gendre 
2000) and with the 1998 June 16.30 UT datapoint of Djorgovski et al. (1998b). The formal 
2o" upper limit to the decay slope of otR < 1.8 indicates that the optical afterglow did not 
exhibit an unusually rapid decay during the first day after the burst. Thus, the hypothesis 
that rapid decay is the cause for the dimness of some optical afterglows (Groot et al. 1998) 
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is not applicable to GRB 980613. 

The spectral index was computed from the NOT R and WHT H data from 1998 June 
13.9 UT, and from the MDM R and I data from 1998 June 14.2. We corrected for a 
Galactic extinction of E(B-V) = 0.09 (Schlegel et al. 1998), i.e., A R = 0.23, Aj = 0.16, 
and A H = 0.05. The resulting spectral indices are /3rj = 0.3 ± 1.1 and /3rh < 1.1. From 
the BeppoSAX X-ray flux reported by Boer & Gendre (2000) and Soflita et al. (2001) we 
determine the spectral index between the R band and the 2-10 keV X-ray band to be 
Aix — 0.59 ± 0.03, 16 hours after the burst, consistent with and the upper limit on /3rh- 
The latter limit translates into a maximum internal extinction in the host galaxy of 0.45 
mag in R in addition to the 0.23 mag of foreground extinction. 

We note that these observations are consistent with a p = 2.2 fireball model (p is 
the power-law index of the electron energy distribution) which in the spherical adiabatic 
expansion (pre-break) regime (Sari, Piran & Narayan 1998) gives (5 = 0.6 and a = 0.9 for 
v < v c . However, given the limited amount of data available, other models are also viable. 

3. The host complex 
3.1. Observations 

We used the Space Telescope Imaging Spectrograph (STIS) aboard the HST to obtain 
CCD images of the environment where GRB 980613 occurred. These images were taken on 
2000 August 20, 799 days after the burst. The total exposure times were 5851 seconds in 
the 50CCD (clear; hereafter referred to as CL) aperture and 5936 seconds in the F28X50LP 
(long pass; hereafter referred to as LP) aperture. The CCD gain was set to 1 e~/ADU and 
the read-out noise was 4.46 e~. We used a four-point STIS-CCD-BOX dithering pattern with 
shifts of 2.5 pixels (= / /127) between exposures. The data were preprocessed through the 
standard STIS pipeline and combined using the DITHER (vl.2) software (Fruchter & Hook 
2002) as implemented in IRAF (v2.11.3)/STSDAS (v2.1.1). The 'pixfrac' parameter was set 
to 0.6 and we selected a final output scale of / /0254 pixel -1 . These observations were taken 
as part of the Cycle 9 program GO-8640: A Survey of the Host Galaxies of Gamma- Ray 
Bursts. The URI for the survey is http://www.ifa.au.dk/~hst/grb_hosts/index.html 
(Holland et al. 2000). In Fig. 2 we show an inverted gray-scale plot of a 5" x 5" excerpt of 
the CL image, centered on the location of GRB 980613. It reveals a chaotic region consisting 
of numerous knots and/or interacting galaxy fragments. 
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3.2. Model fitting 

In order to obtain quantitative measures of the structural and photometric parameters 
of the objects in the host complex we made successive robust fits of a tilted exponential 
disk model convolved by the drizzled combination of a set of 5 x 5 shifted and rebinned 
PSFs calculated by the TinyTim program (Krist 1995) using a subsampling of 5. We fitted 
one object at a time, after having subtracted and/or masked all neighbouring objects. The 
spectral type was chosen to match the CL — LP color of the fitted object as closely as 
possible. A telescope jitter of &!005 was assumed. In order to take proper care of charge 
diffusion in the STIS CCD each of the PSFs in the set was convolved with a 3 x 3 Gaussian 
kernel appropriate for the spectral type used to generate the subsampled PSF. Finally, the 
DITHER software was applied as for the observational data. 

The model fitting is robust in the sense that the parameters are obtained by a maximum 
likelihood (ML) estimation assuming a distribution of the residuals which has broader wings 
than the normal distribution (Press et al. 1992). As an example, the median is an ML 
estimator if the errors are distributed as a symmetric two-sided exponential. We would, 
however, prefer to retain the possibility of having asymmetrically distributed residuals. In 
addition, the deviation from a normal distribution should be small near the peak of the 
distribution. The simplest probability density satisfying these criteria is a convolution of a 
Gaussian and an asymmetric double exponential |7 _1 (7 2 — 5 2 ) exp(— j\x\ +Sx), where 7 > 
and \S\ < 7. The three parameters 7, 5, and the Gaussian a are fixed by use of the maximum 
likelihood estimator on a sky annulus centered on the galaxy to be fitted. The residuals are 
scaled by the square root of the weights produced by the DITHER software. The robust fit 
of the galaxy parameters are obtained by seeking the minimum of minus the logarithm of 
the likelihood function, but this time with 7, 8, and a held at fixed values. 

The minimization is implemented in the IDL programming language by use of the 
downhill simplex method (amoeba) described by Press et al. (1992). In general the errors on a 
parameter determined by a maximum likelihood estimator are not symmetrically distributed 
around the most probable value. The error range corresponding to 2a (if the method of least 
squares had been used) is defined by the parameter values that decrease the maximum of 
ln(likelihood) by 2.0 when the remaining parameters are left free to vary. It was obtained 
by an iterative procedure. In order to obtain convergence for the smallest objects (F, I and 
J in Fig. 3) they were forced to have circular symmetry. 

Fig. 3 is a contour plot of our model fit to the image area shown in Fig. 2. The quality 
of our model is demonstrated in Fig. 4 which shows the residual image smoothed by a 7 x 7 
box. The fit is remarkably good, except for a bar-like structure near the center of the galaxy 
A and some positive residuals in the outskirts of galaxies B and C. In these plots we have 
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extended the labeling introduced by Djorgovski et al. (2000). We use A to refer to the face- 
on disk galaxy (to the north of the GRB, cf. § 3.3). We did not detect the very red object 
D nor the faint object E of Djorgovski et al. (2000) in either the CL or LP images. 

The individual components have apparent diameters ranging from 0'.'04 (0.3 kpc) to 0'/55 
(4.5 kpc). Their small sizes show the need for high-resolution (HST) photometry to resolve 
internal details in GRB host galaxies and their environments. 

3.3. Localization 

We used the discovery NOT i?-band image obtained on 1998 June 13.9 UT to determine 
the location of the optical afterglow in the STIS CL image. The host complex made a 
significant contribution to the flux in the NOT image. Therefore we subtracted the STIS CL 
image after it was transformed, scaled and smoothed to the resolution and intensity scale 
of the ground-based image. Based on four reference sources common to the STIS CL image 
and the ground-based NOT image we find GRB 980613 to be displaced 0'.'20 ± 0'.'05 to the 
west (#o ~ —85°) of a compact, blue galaxy (H in Fig. 3) seen close to the center of the host 
complex, towards a fainter companion (F). Both these objects are embedded in or projected 
onto the large, nearly face-on, low surface-brightness (the cosmological surface-brightness 
dimming amounts to 3.22 mag at z = 1.0969) disk galaxy (A) identified by Djorgovski et al. 
(2000) as the host galaxy based on ground-based images. The high spatial resolution (HST) 
images has allowed us to pinpoint the location of the GRB to the galaxy H which was not 
seen in the ground-based images of Djorgovski et al. (2000). We note that Bloom, Kulkarni 
& Djorgovski (2002a), based on the HST data presented here, also localized the GRB to the 
galaxy H. As shown below, the properties of H are quite different from those of A. 

3.4. Photometry 

We measured the color of each of the objects labeled A, B, C, F, G, H, J, and K in Fig. 3 
from robust model fits to the CL and LP images while keeping the structural parameters 
(scale length, axial ratio, and position angle of major axis) fixed at the values obtained 
from the fit of all model parameters to the CL image. For the color measurements, only the 
positions and magnitudes of the galaxy models were fitted, as it makes little sense to calculate 
the color from magnitudes obtained from different structural parameters (a procedure similar 
to using different apertures). Furthermore, for fixed structural parameters the magnitude 
errors are nearly symmetric so the error in the color can be derived as the square root of 
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the sum of magnitude errors. Note that the CL magnitudes and the structural parameters 
were obtained from unrestricted fits to each galaxy in turn. The rather modest overlaps of 
some of the objects do not justify the complications of a simultaneous multicomponent fit 
(the very robust downhill simplex algorithm is not suitable for problems involving more than 
about 10 parameters). 

The STIS magnitudes CL and LP are calculated on the AB system. We have used the 
zero points 26.386 and 25.291 given by Gardner et al. (2000) for 1 count s -1 in CL and LP, 
respectively. The photometric and morphological parameters are listed in Table 2. For each 
object we present the CL total magnitude, the CL — LP color, the central surface brightness 
A*cl(0) in mag arcsec -2 , the scale length of the exponential disk Rd, the axial ratio b/a, and 
the position angle of the major axis. The errors given in Table 2 correspond to 2a errors, if 
the method of least squares had been used. 

3.5. Star-formation rates 

In the following we will assume that all of the objects in the host complex are at the 
same redshift — that of object A. In this case the separations between the objects indicate 
that the host complex is a compact group. The best fit to the host object H has a disk scale 
length of R d = 0"0360±0 / .'0028 (0.30±0.02 proper kpc). The object is approximately circular 
with an ellipticity of e = 1 — b/a = 0.30. The host object H and its fainter companion F are 
the bluest objects in the field and H has the highest central surface brightness (Table 2). 

As shown in Fig. 1 of Madau, Pozzetti & Dickinson (1998) the power radiated at ~2800 A 
is proportional to the instantaneous Star Formation Rate (SFR) of a model stellar population 
with a stellar birthrate of the form SFR oc exp(— t/r), where r is the duration of the starburst. 
The constant of proportionality is quite insensitive to the details of the past star formation 
history. At z — 1.0969 the rest-frame wavelength of the CL bandpass is very close to 2800 
A. We used the on-line version of the Galaxy Isochrone Synthesis Spectral Evolution Library 
(Bruzual & Chariot 1993; Chariot et al. 1996) and the SYNPHOT package of STSDAS 
under the IRAF environment to calculate both rest-frame and redshifted AB-magnitudes 
for a set of Spectral Energy Distributions (SEDs) for bursts of star formation characterized 
by a Salpeter IMF and an exponentially decreasing SFR with r = 1 Gyr. Specifically, we 
obtained the star- formation rate, the absolute magnitude, M B , the apparent CL magnitude, 
and the CL — LP color for each of 12 starbursts with solar metallicity at ages ranging from 
1 Myr to 10 Gyr. A Galactic extinction of E(B — V) = 0.09 was applied. 

In Fig. 5 we show the specific star-formation rates of the galaxies in the host complex 
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as a function of the observed color, in both the observer's frame (a) and the rest frame 
(b). The two continuous curves represent a cubic spline through the 12 model values. The 
spline curves are used to transform from color to specific SFR and as such are also used 
to infer the errors. Figures 5a and 5b show log(SFR) + 0.4CL and log(S') + 0.4M|> (S is 
the star-formation rate per unit blue luminosity), both as a function of the CL — LP color. 
SFR and S are given in units of .Moyr -1 and A^ yr _1 L^ _1 , respectively, where L* B is the 
rest-frame B-band luminosity of a typical galaxy. Lilly et al. (1995) find a typical magnitude 
of M* = -21.4 (AB) for blue galaxies with z « 1 if {H ,Q m ,Q A ) = (50,1,0). For our 
adopted cosmology of (H ,Q m ,Q A ) = (70,0.3,0.7) this corresponds to M* B = -21.29 (AB). 
The values in Table 2 were calculated using this value. The 2a error bars of CL — LP are 
taken from Table 2. Note that the objects F, J, and K have error bars extending outside the 
full color range of the model SEDs. The end-points of the vertical bars are calculated as the 
curve values of the end-points of the horizontal bars, except for the cases where they extend 
outside the model range. These bars are limited by the extreme values of the burst models. 
The specific SFR with respect to the CL-passband (Fig. 5a) is, as expected, nearly constant 
in the color range 0.0 < CL — LP < 0.5. This is, unfortunately, not the case when the SFR 
is expressed in terms of the rest-frame B-band luminosity (Fig. 5b). 

In the last two columns of Table 2 we present the specific star-formation rates in units 
of MqYT^ 1 L*^ 1 for eight of the objects in the host complex, as well as M B — M B for three of 
these objects, including the host object. The intrinsic luminosities for the rest of the objects 
are not given due to their very large uncertainties. The star-formation rates listed in Table 2 
and the curves shown in Fig. 5 assume that there is little or no intrinsic extinction in the 
galaxies of the host complex. The star-formation rates determined from the optical images 
therefore represent lower limits to the true values. 

The star- formation rates per unit blue luminosity in the compact host object (H) and 
its companion (F) are the highest of the objects listed in Table 2. If H is a galaxy in its own 
right then it is sub-luminous with Lb = (0.05 ± 0.01)L^. However, this value is somewhat 
uncertain since M B is highly correlated with both the slope of the faint end of the galaxy 
luminosity function and its normalization (Lilly et al. 1995). The specific star-formation 
rate S of the compact host is ~ 17 A4 Q yr^ 1 L* B ~ X , the highest value of any known GRB 
host galaxy, as estimated from the 2800 A flux. The star-formation rate density inside the 
disk scale length is xs 12 M.® yr _1 kpc~ 2 , consistent with values (1-1000 A4® yr -1 kpc~ 2 ) 
characterizing a star burst (Kennicutt 1998). We stress that this is lower limit to the true 
star- format ion rate because of the unknown UV extinction in the host galaxy (see also Chary 
et al. 2002). 
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4. Discussion 

The dim optical afterglow of GRB 980613 would have gone undetected under slightly 
less favorable search and observing conditions and, in this case, would have been classified as 
a 'dark burst' (Fynbo et al. 2001; Lazzati, Covino & Ghissellini 2002; Reichart & Yost 2001). 
Part of the reason that the optical afterglow was faint may have been dust extinction in the 
host galaxy. However, we have presented evidence that the afterglow suffered only at most 
0.45 mag of extinction along the line of sight (at a redshift of 1.0969 the rest-frame reddening 
in the host is Ay < 0.25) in the host galaxy plus 0.23 mag of foreground extinction. 

Another explanation for the faintness of the optical afterglow may be that the burst 
was collimated but with the line of sight outside the collimated beam. This results in a 
much fainter peak afterglow magnitude while the late time light curve behaviour is similar 
to that observed if the line of sight is inside the collimated beam (Huang, Dai & Lu 2000). 
As observed, the burst may in that case be expected to be weaker in 7-rays than in X-rays 
(Nakamura 2000). 

The positional coincidence between GRB 980613 and the faint compact blue object (H) 
with a high surface brightness and specific star-formation rate provides strong evidence in 
favor of the association of GRBs with star formation. Not only did GRB 980613 occur in a 
chaotic and possibly interacting region (Djorgovski et al. 2000), the STIS images presented 
here also reveal that the GRB occurred in what is probably the most active site of star 
formation in the entire complex. 

It is interesting to note that GRB 980613 did not occur in the region of highest inferred 
total star-formation rate (object A) but in the region of the highest specific star-formation 
rate (objects H and F) (although the total star- formation rate of H is also relatively high). 
One may speculate that GRBs trace regions of high specific optical star-formation rate, 
either because the total star-formation rate is in fact dramatically higher than that inferred 
from the UV flux due to dust (as in GRB 980703, cf. Berger, Kulkarni & Frail (2001)) or 
because the IMF could be biased towards high-mass stars in regions experiencing its first 
intensive starburst, possibly leading to a much higher rate of massive-star formation. We 
caution however that such conclusions should not be based on one system and note that 
GRB 990705 (Andersen et al. 2002) occurred in a region of very little star-formation activity 
in an evolved galaxy. 

JH and AO J thank Rolf Stabell and Sjur Refsdal for their hospitality at the Center for 
Advanced Study in Oslo where the identification of the optical afterglow of GRB 980613 
was made. We thank Enrico Costa and Eliana Palazzi for early information that facilitated 
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Fig. 1.— i?-band images obtained with the NOT on 1998 June 13.9 UT (left) and 17.9 UT 
(right) which led to the identification of the optical afterglow of GRB 980613 marked by an 
arrow. The images are V x 1'. North is up and east is to the left. 
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Fig. 2.— Inverted gray-scale plot of a 5" x 5" excerpt of the HST/STIS CL (50CCD clear 
passband) image of the host galaxy of GRB 980613 and its immediate surroundings. The 
location of the optical afterglow of GRB 980613 is marked by a cross, the 2a error ellipse 
(dashed), and an arrow. 
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Fig. 3. — Contour plot of a model fit to the 5" x 5" area shown in Fig. 2. The location of the 
optical afterglow of GRB 980613 is marked as the 2a error ellipse (dashed) and is marked 
with an arrow. A blue compact object H is evident in the center of the image. The objects 
A, B, and C were identified in ground-based images (Djorgovski et al. 2000). We identify F, 
G, H, I, J, and K as well. 
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Fig. 4. — Residual image obtained by subtracting the model fit shown in Fig. 3 from the 
HST/STIS CL image. It has been smoothed by a 7 x 7 box in order to suppress the small 
scale noise. The area shown is exactly the same as the one shown in Fig. 3. The crosses 
indicate the positions of the model galaxies, and the labeling is the same as in Fig. 3. The 
brightness varies from —9% to +12% of the average sky background. The fit is remarkably 
good, except for a bar-like structure near the center of galaxy A and some residuals close to 
galaxies B and C. 
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Fig. 5. — Specific star-formation rates of the objects in the host complex expressed in the 
observer's frame (a) and in the rest-frame (b), respectively. Panel (a) shows log(SFR)+0.4CL 
as a function of the apparent color CL — LP (CL is the 50CCD clear passband, LP is the 
F28X50LP long passband). As expected, the SFR scales with the apparent magnitude CL. 
Panel (b) shows log(S') + 0.4M^ as a function of CL — LP. Note that S scales with the 
chosen reference luminosity L* B . The curves plotted represent a model with an exponentially 
decreasing SFR with r = 1 Gyr and a range of ages between 1 Myr and 10 Gyr. SFR 
and S are in units of A^ yr _1 and .Moyr -1 !,^ -1 , respectively. The theoretical curves in 
(a) and (b) are used to translate the measured CL — LP colors with error bars into specific 
star- formation rates with error bars (it is not a fit). The 2a error bars of CL — LP are taken 
from Table 2. 



Table 1. Log of observations and photometry of the afterglow of GRB 980613 



Date 


Filter 


Exposure time 


Telescope a 


FWHM 


Brightness 13 


Afterglow 


(June 1998 UT) 




(sec) 




(arcsec) 


(mag) 


(mag) 


13.89 


R 


600 


NOT 


0.7 


22.81 ±0.15 


23.19 ±0.20 


13.89 


K> 


5 x 450 


CAHA 


1.2 


18.8 ±0.3 




13.90 


H 


2400 


WHT 


1.2 


> 20.9 




14.20 


I 


6 x 500 


MDM 


0.96 


22.64 ±0.08 


23.13 ±0.20 


14.24 


R 


4 x 500 


MDM 


1.14 


22.96 ±0.09 


23.53 ±0.20 


14.90 


K> 


4 x 600 


CAHA 


1.2 


> 18.5 




15.19 


I 


4 x 900 


MDM 


0.99 


23.25 ±0.15 


24.1 ±0.4 


17.91 


R 


600 


NOT 


0.68 


> 24.4 





a NOT is the 2.56-m Nordic Optical Telescope on La Palma, equipped with ALFOSC 
(1998 June 13.89) and HiRAC (1998 June 17.89); CAHA is the 3.5-m telescope at Calar 
Alto, equipped with OMEGA-Cass; WHT is the 4.2-m Willian Herschel Telescope on La 
Palma, equipped with CIRSI; MDM is the 2.4-m Hiltner telescope of the MDM Observatory 
on Kitt Peak, equipped with a thinned, back-side illuminated SITe CCD. 

b The R and I data were calibrated to the Cousins system using Landolt (1992) standard 
stars (Halpern & Fesen 1998). The quoted magnitudes refer to a 1" radius aperture. 

c The R and / magnitudes after subtracting scaled and transformed HST/STIS images to 
correct for the contribution from the host environment. 
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Table 2. Photometry, morphology and star-formation rates of the host complex 



Object 


CL 


CL - LP 


MCl(0) 




b/a 


P.A. 


log(S) 


M B - Mg 




(mag) 


(mag) 


(mag arcscc~ 2 ) 


(mas) 




(°) 


(Afoyr- 1 !^- 1 ) 




A 


24.82 +° ™ 


0.33 ±0.10 


24.02 


277 +f 4 


0.80 


6 


0.82^? 


1.3 ±0.3 


B 


25.80 t°;°S 


0.48 ±0.10 


22.62 


91.9 ±1* 


0.39 


104 


0.64t° il 




C 


9 C OO +0.09 


0.62 ±0.08 


22.43 


83.0 tfl 


0.34 


160 


o.i5tr 6 6 9 




F 


27.34 +°;|° 


-0.07 ±0.30 


22.39 


41.0 in 


1.00 








G 


27.11 ±°;it 


0.33 ±0.19 


21.44 


29.3 ±1* 


0.31 


52 


0-82±g;2| 


3.6 ±0.5 


H 


26.08 IHo 


0.09 ± 0.08 


20.86 


36.0 til 


0.60 


122 


1 r.q + 0.21 


3.3 ±0.3 


I 


28.10 t° f 7 




21.55 


19.5 til 


1.00 






J 


28.06 + 02 7 


0.65 ±0.41 


22.28 


27.9 


1.00 




-0 1+ ' 9 
u - i -1.6 




K 


26.03 ±°iS 


0.57 ±0.23 


24.40 


188 tf 5 


0.07 


16 


4+ ' 4 
u - 4 -2.0 





Note. - 



The meaning of the quoted errors is discussed in §§ 3.2 and 3.4 



